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abstract 

In this paper, we consider the discretization 

y{t + e) = y{t-e) + 2e{l-v{tf), 

e > a small parameter, of the logistic differential equation y' — 1 — y"^, 
which can also be seen as a discretization of the system 




This system has two saddle points at A ~ (1,1), B = (—1,-1) and there 
exist stable and unstable manifolds. We will show that the stable manifold 
of the point A = (1,1) and the unstable manifold W~ of the point 
B — (—1,-1) for the discretization do not coincide. The vertical distance 
between these two manifolds is exponentially small but not zero, in particu- 
lar we give an asymptotic estimate of this distance. For this purpose we will 
use a method adapted from the paper of Schafke-Volkmer ^ using formal 
series and accurate estimates of the coefficients. 

Keywords : Difference equation ; Manifolds ; Linear operator ; Formal so- 
lution ; Gevrey asymptotic ; Quasi-solution 



1 Introduction 



We consider the logistic equation 

y' = l-yitf, (1.1) 

whose solutions are y{t) = tanh(t + c). The discretization of this equation by 
Nystrom's method, which consists in replacing the derivative by the symmetrical 
difference, gives the recurrence : 

yn+i = yn-i + 2e{l-yl). (1.2) 
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Fig. 1 - This figure represents y„ as a function of ne^. 



With the initial conditions yo = 0, yi = e, we calculate the discrete solution y„ 
with 1600 iterates and e = ^ (see FIG.l). We observe that the discrete solution 
rejoins quickly the level y = 1, it stays close to this point for a relatively long time 
before it leaves its neighborhood, then for the little while the even numbered points 
follow one curve and the odd numbered follow another curve, the two curves meet 
at the level y = — 1, the discrete solution stays close to this point for a relatively 
long time, then it remakes the same cycle. 

Now letting Un = y2n, Vn = y2n+i, wc obtain a recurrence of first order in the 
plane 

{Un+l,Vn+l) = ^{Un,Vn) (1.3) 

where the diffeomorphism <I> : IR^ i-^ IR^ is defined by 

ui = u + 2e{l -v^), 
vi = V + 2e{l - uj). 

This is a discretization of the following system of differential equations 

(1.4) 

We notice easily that the set (B = {{u, v)\{u — v){u^ + + f ^ — 3) = O} is an 
invariant set for this system. The system (11.41) has two saddle points in ^ = (1, 1) 




and B = (—1,-1) and corresponding stable and unstable manifolds lie on its 
invariant set ^. 

The stable manifold at A and unstable manifolds at B are part of set { {u, v)\u = 
u}, where as the unstable manifold at A and stable manifolds at B are part of set 
[{u,v)\u'^ + uv + v'^ = 3^. The stable manifold of A coincides with the unstable 
manifold of B (See FIG.2). 

For the discretized equation (|1.3p . these manifolds still exist P], let , W~ 
denote the stable manifold at A and the unstable manifold at B respectively, 
and W~ unstable manifold at A and stable manifold at B respectively, and 
Wj^ do not coincide an more (See FIG. 3) as we want to show. 

In the paper [2], after introducing the notion of length of the first level li{e) = 
2eni(e), where ni(e) = inf {n G lN\y2n(e) +y2n+i(£) < O}, Fruchard-Schafke had 
shown that there exist a postive constant K such that 

A(e) < exp - ^^^^) , as e \ 0, 

^i(^) > ^ + 0(1)' 

where A(e) denote the distance between the sets H S and W~ n 5 in the sense 
of Hausdorff and S = [{u, v)\ — l<n + u<l}. 

With the initial condition yo = 0, yi = e, the length of the first level satisfies 
h{e) = — ^(l + o(l)) log(A(e)) as e \ Of^j- They also showed that there are two 
families of entire functions yf : C i — > C, solutions of the difference equation 

y{t + e)=y{t-e) + 2e{l-y{tf), (1.5) 

and the functions t i— > [yf{t), yf{t + e)) , t ^ (^y^{t),y^{t + e)) provide parame- 
trization W^it) of for t G [— 1, oo[ respectively w~ (t) of W~ for t g] — oo, 1]. 
In this work we will prove 

Theorem 1.1. There exist a constant a with 1.2641497 < a < 1.2641509 and 
eo > such that for < e < Eq 

, , , , , 47racos(-t + vr) ^2 f 1 ^2 \ 

where Dists is defined as the vertical distance between the stable and unstable 
manifolds. 



In order to show the theorem 11.11 we start with the construction of a formal 
power series solution in d whose coefficients are polynomials in u = tanh((it/e), 
afterwards we will give asymptotic approximations of these coefficients using ap- 
propriate norms on the spaces of polynomials. The next step is to construct a 
quasi-solution i.e. a function that satisfies the equation (II. 5p except for an expo- 
nentially small error, then we show that this quasi-solution and the exact solution 
of equation (II. 5p are exponentially close. 
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Fig. 2 - The stable and unstable manifolds for the logistic differential equation. 

Finally we give an asymptotic estimate of the distance between the stable 
manifold of A and the unstable manifold of B and we show that this distance is 
exponentially small but not zero, thus completing the proof of the theorem 11.11 




Fig. 3 - The stable and unstable manifolds for the difference equation. 



2 Formal solutions 

We consider the difference equation 

yit + e) = y{t-e) + 2e{l-yitf), (2.1) 

where e > is the discretization step and y{t) — > 1 as t — > +oo. Our first aim is 
to transform (12. 1|) in such a way that the new equation admits a formal solutions 
whose coefficients are polynomials. We define u and A{d, u) such that 

u := tanh(-t), (2.2) 
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ye{t) =A{d,u), (2.3) 

where d := Yln=L '^n^^ is a series to be determined. When these variables are 
substituted in (12. the following equation satisfied by A{d,u) is obtained : 



A[d, T+{d, u)) - A{d, T-{d, u)) = f{e, Aid, u)) , (2.4) 



where 



f{e,x):=2e{l-x'), (2.5) 
„4- / , ^ u + tanh((i) 

r+ := r+ d, u) = 2.6 

^ ^ l + ntanh(d)' ^ ' 

For small e one can construct a formal expansion in powers of of the form 
^2n+i(^)c^^", where A2n+i are polynomials which satisfy Ai{u) = u and 
^2n+i(±l) =0, forn> 1. 

The existence of a such formal solution is only possible if d and e are coupled 
in a very special way. Indeed, suppose that there exists a such formal power series 
solution. We differentiate (|2.4I) with respect to u and obtain at u = 1 



Because A2n+i(l) = for n > 1 implies A{d,l) = 1. Thus with ^^(e, 1) = -4£ 



dx 

and ^(d, 1) = 1 + 0((i) 7^ 0, the above equation implies 



tanh(<i) _l^.^j^(2d). (2.8) 



1 - tanh((i)2 



Theorem 2.1. Lei e and d be coupled by £ = ^ sinh(2d). Then {2.4^ has a formal 
power series solution that can be written in the form 



A[d, u)=u + Y^ ^2n+i(n)(i2", (2.9) 

n=l 

where A2n+\{u) are odd polynomial and A2„+i(l) = ^2n+i(— 1) = for all n. 

Proof. We will use the Induction Principle for even n to show that there exist 
unique odd polynomials Ai, A^, A^...An+i such that 

n 

Z{d,u) = Ak+id\ (2.10) 

fc=0 

k even 



5 



satisfy 

Z{d,T+{d,u)) - Z{d,T~{d,u)) =f{e,Z{d,u)) modd''+^. (2.11) 
For n = 0, we put Ai{u) = u and Z{d, u) = u and obtain 
Z{d,T+{d,u)) - Z{d,T-{d,u)) =T+{d,u) -T-{d,u) = {2 - 2u'^)d + 0{d^). 
and 

f{e,Z{d,u)) =2e{l-u^) = {2 - 2u^)d + 0{d^). 

This gives 

Z{d,T+{d,u)) - Z{d,T-{d,u)) =f{e,Z{d,u)) mod d^ . 

Now suppose that for some even n aheady Ai, A^, Ar,...An+i have been found 
with the above properties. We have to construct An+z- First, we show that Z{d,u) 
satisfies (|2.11l) even modulo d""*"^. To this purpose, let 

Z{d, r+) - Z{d, T-) = f{e, Z{d, u)) + i?„+2(n)d"+2 + 0(d"+3) (2.12) 
we replace d by —d. Using that Z is even and e is odd in d, we obtain 

Z{d,T-) - Z{d,T+) = -2e(l - Z{d,uf) + i2„+2(n)d"+2 + 0(d"+3), 
this gives 

Z(d,T+) -Z(d,r~) =2e(d)(l-Z(d,7x)2) -i?„+2(n)((i)"+2 + 0(d"+3). 
With (|2.12p this implies Rn+2iu) = and consequently 

Z{d,T+) -Z{d,T-) = f{e,Z{d,u))+Rn+z{u)d''+^ + 0{d^+^). (2.13) 
We want to construct An+?,{u) such that 

Z{d, T+) - Z{d, T-) = f{e, Zid, u)) + Oid^+^), (2.14) 

if we put 

Z{d, u) = Z{d, u) + ^„+3(u)d"+2^ (2.15) 
To this purpose we use again Taylor expansion 

f{e,Z{d,u)) = f{e,Z{d,u)) -4uAn+3{u)dJ'+^+ 0{d^+^), 

Z{d,T+{d,u)) = Z{d,u)+An+3iu)d''+^ + {l-u^)^^^{u)d''+^ + 0{d^+^), 

Z{d,T-{d,u)) = Z(d,n)+A„+3(^^)d"+'-(l-n2)^^^(^.)(i"+3 + (9(d"+4). 
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With (I2l3]) . this gives 

Z{d,T+) -Z{d,T-) = f{e,Z{d,u)) + Rn+3{u)d^+^ + 0{d^+^), (2.16) 

where 

Rn+3{u) = 2(1 - u^)^^^{u) + AuAn+3{u) + i?„+3(n). 

We see that (|2.14l) is satisfied if only if 

dA 

2{l-u^)^^{u)+4uAn+3{u) + Rn+3{u)=0. (2.17) 
This equation has a unique odd solution which is given by 

A„^3iu) = -{l-u')£^^±^dt. 

We can prove that this solution is polynomial. For that, it is necessary that 
Rn+3{u) and R'^_^_-^{u) vanish into 1 and -1. Indeed, if we take (I2.13P with u = 1, 
we obtain 

i?„+3(l)d"+3 = Z{d, T+{d, 1)) -Z{d, T-{d, 1)) -f{e, Z{d, 1)) + 0(d"+^). (2.18) 

Since, T+{d, 1) = T-{d, 1) = Z{d, 1) = 1 and f{e, l) = 0, we obtain 

i?„+3(l)d"+' = 0{d^+^), (2.19) 

therefore -R„+3(l) = 0. In order to show that R'^^^{1) = 0, we derive formally (|2.13p 
and replace u by 1. Using T+{d, 1) = T-{d, 1) = Z{d, 1) = 1 and §^(e, l) = -4e, 
we obtain 

<,+3(l)d"+3 = 1) [-^{d, 1) - -^{d, 1) + 4eJ + 0(d"+4). (2.20) 

By choice of d, cf Thus we have ^{d, 1) - ^{d, 1) = -2sinh(2d). with 

(2,8) we obtain 

<+3(l)d"+3 = (2.21) 
therefore -^^^3(1) = 0. As formula (|2.13|1 shows that RnJ^^iu) is odd, then i?„_|_3(— 1) = 

i?;+3(-l) = 0. 

The first polynomials An+i{u) with even n are given by 



n 





2 


4 


6 


An+l{u) 


u 


u — V? 


(l-n2)(|tx-fn3) 









We introduce the operators C2,C,52,5 defined by 
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C{Z){d,u) = ^{Z{d,T+l) + Z{d,T-'2)), 
S{Z){d,u) = \{Z{d,T+-2) - Z{d,T-"2)), 

(2.22) 

C2{Z){d,u) = ^{Z{d,T+) + Z{d,T-)), 

S2{Z){d,u) = \{Z{d,T+) - Z{d,T-)). 

where T~^2 = r+(|,ii), T~2 = T~{^,u) and Z{d,u) is a formal power series of d 
whose coefficients are polynomials. We rewrite equation (|2.4|) as 

S2{A){d,u)=e{l-A{d,u)^). (2.23) 
3 Norms for polynomials and basis 

In the sequel we denote : 

- V the set of all polynomial whose coefficents are complex, 

- Vn the spaces of all polynomials of degree less than or equal to n, 

- Q:= {Q{d,u) = En=oQn{u)d'', where Q„(n) G Vn, for all n G N} . 

Proposition 3.1. If we define the sequence of the polynomial functions r„(u) 
by tq{u) = l,ri(n) = u,Tn+i{u) = ^DTn{u), where the operator D is defined by 
D := (1 — u'^)-§^, we have 

1. T+{d,u)=Y.n=^rn+l{u)d^, 

2. Tn{u) has exactly degree n, 

3. r„(tanh(z)) = ^(A)"-^(tanh(z)). 

Definition 3.2. Let p G P„. As ro(n), ri(u), r„(u) form a basis ofVn, we can 
write p as 

n 

p = ^afcrfc(u). 

fc=0 

Then we define the norm 

n 
i=0 



Theorem 3.3. Let n,m be positive integers and p G Vn, q G Vm- The norms 
II ||„ of the above definition have the following property : 
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1. \\Dp\\n+l < n\\p\\n- 

2. Ifp odd we have \\p\\n < ||-Dp||„+i. 

3. There exists a constant Mi such that \\pq\\n+m < -^illpllni'zllm- 

4. There is a constant M2 such that \p{u)\ < M2 (f)'" \\p\\n {-I < u < 1). 

5. There is a constant M3 such that for all n > 1 with p{—l) = p(— 1) = 0, 

p 



T2 



<Ms\ 



4 Operators 

In this section we will use definitions and results adapted from [5] by replacing 
f by TT. 

Definition 4.1. Let f be formal power series of z whose coefficients are complex. 
We define a linear operator f{dD) on Q by 



00 n 



/(dZ))Q = (4.1) 

n=0 i=0 

Where f{z) = TZ^ h^" (^nd Q = En=oQniu)d^ G Q- 

By the above definition and (1) of the proposition 13.11 we can show that 

Q{d,T+{6d,u)) = {eicp{6dD)Q){d,u), for Q G Q and all 6 £ C. (4.2) 

With K22\\ this implies 

C{Q){d,u) =cosh(fL>)g, S{Q){d,u) = sinh(f L')Q, 
C2{Q)id,u) = cosh{dD)Q, S2{Q){d,u) = smh{dD)Q, 

for polynomial series Q in Q. We denote ||Q||^ = ||Qn||„- 

Theorem 4.2. . Let f{z) be formal power series having a radius of convergence 
greater than vr and let k be a positive integer. There is a constant K such that : If 
Q is a polynomial series having the following property 



IIQIIn < 



for n < k 

M{n - k) !7r~" for n > k 



where M independent of n and Q £ Q then the polynomial series f{dD)Q satisfies 
\\f{dD)Q\\n < 



for n < k 

MK{n - k) Ivr"" for n > k 
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We define the operator C ^ by C ^ = g{dD), where g{< 



cosh(2;/2) 

Theorem 4.3. There exists a positive constant K such that, if Q is a polynomial 
series such that Qn odd, \\Q\\n = for n < k for some positive integer k and 

||dL»Q||„ < M(n- A;)!7r-" forn>k 

where M independent of n and Q E Q, the polynomial series C~^{Q) satisfies 



\C-\Q)\l<MK7T- 



n ! for k = 1 

(n-l)!log(n) fork = 2 
^ (n - 1) ! fork>3 



Theorem 4.4. We consider a polynomial series 



Qa{d,u) = an{n-l)\y-j Tn{u)d 



n=l 
n odd 



where an = 0{n ^) as n ^ oo with some integer k >2. Let a := YlnLi 
the coefficients {C~^((5a)}„ of C~^{Qa) satisfy 

{C-\Qa)}^ - a(n - 1) !(^)"~Vn = 0((n - A:) Itt"-) , 



len 



as n ^ oo for odd n. 

Theorem 4.5. Let k,l,p,q be positive integer with p > k and q > I. Define m as 
the minimum of k + q and I + p. Then there is a constant K with the following 
property : 

If P and Q are polynomial series such that \\P\\n = for n < p, \\Q\\n = for 
n < q and 

||P||„ < Mi(n-A;)!7r-'*, forn>k, 
IIQIIn < M2(n-0!7r-", forn>l. 

Then 



\\PQ\\n < KMiM2{n-m)\Tr-"', for n > p + q. 



Theorem 4.6. Let Qi{d,u) be a convergent polynomial series which is even with 
respect to both variables and has constant term 1. 

Let Q2{d,u) = d!^{l — u'^)Qi{d,u) and P{d,u) = S{Q2)/C{Q2)- Consider the linear 
operator defined by 

C{Q) = S{Q) - P{d, u) ■ C{Q), Q€Q. (4.3) 
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Then, there exist a constant K with the following property. If Q is an odd polyno- 
mial series with odd coefficients Qn{u) satisfying (5n(l) = for all n, \\C{Q)\\n = 
for n < 6 and 

\\C{Q)\\n < M{n - 6) Ivr-" for n > 6, 

then also 

WdDQWn < KM{n - 6) Ivr"" for n > 6. 



5 Asymptotic approximation of the coefficients of the 
formal solution 

The objective in this section is to construct an asymptotic aproximation of the 
coefficents of the formal solution (|2.9p . It will turn out to be convenient to consider 
the new series B{d,u) = eA{d,u), 

oo 

B{d,u) = eA{d,u) = ^ An+i{u)ed"- -\- eu, 

n=2 
n even 

this gives 

oo 

B{d,u)= Bn{u)d^ + eu, 

n=3 
n odd 

where Bn{u) are odd polynomials. Furthermore we have Bn{u) € Vn and Bn{l) = 
Bn{—1) = for all n. The new equation for B{d,u) is 

S2iB){d,u) =e'^ -B{d,uf. (5.1) 

We saw in section 2 that the series A was a formal solution of the starting equation. 
For the moment, nothing is known about the norms of its coefficients, but we will 
show that this series is Gevrey-1, more precisely ||^||n = 0(n!7r~"). This enables 
us thereafter to construct a quasi-solution. To this purpose, we will prepare the 
equation (|5.1|) so that we can construct a recurrence. 

We start with the decompostion of series B{d,u) in the form : 

B = U + F, where U = eu + {u - u^)d^ + (— - —u^ + 2u) d^ . (5.2) 

V 3 3 / 



Next, we define 



G: = QF, 

J-.^fl. (5.3, 
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where 



\ 2 6 / 

/47 o 31 58 4 104 2 -,\ ,s 
V 6 6 15 45 / 

Qi(d,n) = (l-u2)d2 + ^(n2-u^)d^ (5.4) 

The choice of Q and Q\ and by using the properties of the operators 5, C, we 
will be able to rewrite the equation (|5.ip in the form : 



eo(d,n)£(J) = ei(d, u)C( J) + e2(d, n)F + e3(d, n)^^ 
+e4(d, u)C2{F) + e5{d, u) 

where C is the operator defined in (14.3(1 and ei{d,u),i = 0, ...,5 are convenient 
known convergent series in d and u which will be thereafter given. 

The left hand side of this equation is series with leading term 1 multiplied 
by the invertible operator C applied to the series J. The right hand side is an 
expression of F and J multiplied by known convergent series ei{d,u),i = 1, ...,5. 

U, Q, Qi were chosen so that the series ei{d,u),i = 0..5 are of a rather large 
order in d, this makes the second term smaller than the right hand side. This 
property will be useful to construct a recurrence on n and to reverse then the 
operators C,C, which makes possible to estimate the coefficients of series C{G) and 
finally the coefficients of the formal solution of the equation (|5.ip . 

More precisely, we insert (|5.2|) into the equation (|5.ip and find 



S2{F) = -F^ -2U F + -U^ -S2{U). (5.5) 

We define X by 

X ■= V{d, u) ■ SC{G) + W{d, u) ■ C'^(G) (5.6) 

where 

V{d,u): = l + 2^2^2_(7 4 _ 3^2)^4 _ ^^^6 _ 221^4 +2|^2_ 41X^6^ 

o V o y 15 15/ 

W{d,u): = 2eu. 

Then, because of G = QF, we have X{G) = V ■ SC{F Q) +W -C^iFQ). Using 
the formulas 

C2 = 2C2 - Id, 
S2 = 2SC, 

(5.7) 

S2{FQ) = S2iF)C2{Q)+S2iQ)C2iF), 
C2{FQ)=C2{F)C2{Q)+S2{Q)S2{F). 
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We obtain 

X = W-S2{FQ) + ^W-(C2{FQ) + Qf'^ 

= 1{VC2{Q) + WS2{Q))s2iF) + l(y 52(Q) + W C2iQ))c2iF) 
+ \WQ-F 
With (|5.5p this imphes 

X = ]^WiF + ]^W2C2{F) - ]^{VC2{Q) + WS2{Q))F^ + i2{e,u), (5.8) 

where 

Wi = -2U{VC2iQ) + WS2iQ)) +WQ, 
W2 = VS2iQ) + WC2iQ), 

^l{e,u) = ^iVC2{Q) + WS2{Q))ie^ - U^-S2{U)). 

The lowest power of d in Wi and W2 is the eleventh Wi = 0{(f^) = W2, fi{e, u) 
is analytic. On the other hand C{G) = Qi J, this implies 

X{G) = V ■SiQiJ) + W -CiQiJ). (5.9) 

If we use the product formulas 

S{QiJ)=S{J)C{Qi) + S{Qi)C{J), 

(5.10) 

C(Qi J) =C(J)C(Qi)+5(Qi)5(J). 

We have 

X{G) = (vCiQi) + WS{Qi)y{J) + (vS{Qi) + ^C(Qi))c(J). (5.11) 
With dSSl) and (^3]) . this imphes 

ViC{J) = -{Ws + PVi)C{J) + UviF - ]^(yC2{Q) + WS2{Q))F^ 

+ \w2C2{F) + ^l{e,u), (5.12) 

Vi = VC{Qi) + WS{Qi), 

= VS{Qi) + WC{Qi), 

P = S[{l-u^)(f)/C{{l-u^)(f) = -ud + 0{(f), 

C{J) = S{J)-P-C{J). (5.13) 
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We divide (|5.12p by (1 - n^)(i^, this implies 



(l-ii2)d2 ^ ' (l-n2)(i2 ' ' 2(1-^2)^2 2(l-n2)d2 

■ + (5.14) 



2(l-n2)(i2 ^ (l--u2)d2' 
where F and J are coupled by the equation QiJ = C{Q F) = C{G). 
Theorem 5.1. With the above notations, we have the following estimation 

\\dDC{G)\\^^^ = 0(^{n - 7) !^-") as n ^ oo. 



Proof. We set 

_ \\dDC{G)\\^^^ 

(n-6)!7r-" " ^ ' 

We must show that = 0{n~^). In the sequel we will use the convention : If 
an,n = 0, 1, .... is any sequence of positive real numbers, then 

:= max(ao, oi, ...a„), for all n > 0. (5.16) 

In this proof Ki, K2, ■■■ will always denote constants independent of n. Theorem 
14.31 gives 

||G||„<irie+(n-l)!7r-", (5.17) 
llG^II < i^2/n(ra- 10)!7r~", forn>18, (5.18) 



where 

n-9 



fn = Tetet ^ ^'-^y-^""-'- forn>18. 

i=9 



(5.19) 



(n - 10) ! 

Using theorem 13.31 (5), we find 

||F|L<i^4e++i(n + l)!7r- 

< i^3/n+4(n - 6) Ivr-. 

Using theorem 14.21 we obtain 

||C2(F)||^<K5e++i(n + l)!7r-«. (5.20) 

Observe that the convergent polynomials series Wi/2{l—u'^)(P , W2/2{l—u'^)d? 
begins with (f and [V C2iQ) + W S2iQ)) /2{1 - u'^)(f begin with 1. Using theorem 
14.51 we thus obtain 

^(^^C.(F)|[<AV:_,(„-7)!.-», (5.21) 
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Wi 



2(1 - n2)d2 
[VC2{Q) + WS2{Q)) 



<K,e+_,{n-7) Ivr"", 



2{l-v?)d? 



(5.22) 
(5.23) 



On the other hand 

||QiJ||„ = ||C(G)||„<e„(n-6) Itt'^ 
We apply the theorems 13.31 (5) and 14.21 and find 

||C(J)||„ <i^9e„+2(n-4) Ivr-'^. 
Because the convergent polynomial series (W3 + PVi)/(l — u'^)d? = 0{d^), then 
{W3 + PV1) 



(1 - u'^)(P 



-C{J) 



<Kioe+„i(n-9) Ivr-". 



(5.24) 



With a crude estimate of the convergent terms ^[e, u) the inequalities (|5.2ip - (|5.24p 

gives 



2{l -u^)d? 



C{J) < Kn{l + e+_i + /++4)(n - 6) In 



The multiplication by the convergent terme 2(1 — u'^)(P / Vi only changes the 
constant Kn. Since theorem 4.5 applies to the operator C defined in (14.31) . we 
obtain 



\\dDj\l < K,2{1 + e+_i + /+.4)(n - 6) !^-". 
This with theorem 14.51 gives via 

dDC{G) = dD{QiJ) = J • (dDQi) + Qi ■ {dDJ), 

the estimate 

\\dDC{G)\l < K,,{1 + e+_3 + - 8) l^"". 

By equation (|5.15ll . we obtain 

en-i < -(l + e+ 3 + /, 



(5.25) 



(5.26) 



n 



n-3 "T Jn+2J- 



(5.27) 



Lemma 5.2. Under the condition h5.21\) . we have e„ = 0{n ^) as n ^ 00. 
Proof. Let Ki >9 lefg an arbitrary number. We assume that 

Ki{n + p) ! 



en < 



(n - 1) \ip + 10)1' 



for 10 < n < - 3, 



(5.28) 
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with J9 > -1, iV > 14. This gives for 13 < n < TV 

in - 8) < 2e, K,^^-^-^ + g ^— . 

Using the inequahty, 

n-8 

^(i+p)!(n + p-i + 2)! < {p+W)\{n + p-7)\, 

1=10 

we obtain 

^"+^-^^ (p+10)!(n-8)! -^ ^p+10)!(n-4)! 13 < n < iV, 

with a constant K2 depends on ki, independent of p. The asumption (|5.27p of the 
lemma yields 

< + K, + K.,)^Jf±I^^ f„. 13 < „ < N, 



this implies 

Cn < —1 ^ — TT, for 12 < n < iV - 1, 

n + l(p + 10)!(n-l)! " " 

with a constant K3 only depending upon Ki. Now we choose Nq > 12 so large 
that ^ < Ki and then p so large that (|5.28|) holds. Then 

^ ^im^/"^' iM for 12 < n < TV - 1, with N > No- 

n + 1 + 10) !(n — 1) ! 

Since Ki >9 !eio is an arbitrary number, we have shown for any p > —1 , taht 

„/(n+p)!\ 

implies that 

/(n + p-l)!\ 
"" = ^V (n-1)! J 

Consequently 

Cn = 0[n~^), as n — > 00 

Thus we have shown that 

\\dDC{G)\\^^^ = 0{{n - 7) Ivr-") as n ^ 00 . 
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Let E := C{G) . Like G, the polynomial series E is odd in d and the coefficients 
are odd in u. We partition them 



)n— 1 f i \ 

Tn{u) + l3n-2{n-2,)\\^-j Tn-2{u) 



i \ n— 5 

+ 7„_4(n - 5) ! ( - r„_4(n) + £;„_6, 



vr 



(5.29) 



for odd n > 7, where a„ and /3„ are real number and also E have at most degree 
n for all n. For the whole series E this is equivalent to 



C{G) = E = Ei+ (fE2 + d^^4 + d^^, 



(5.30) 



where 



+00 



E^ 



Eo 



E 



^an{n-l)\i-\ Tn(u)d 



n=7 

+00 



n-1 



5^/?n(n-l)! 

n=5 
n=3 

^:B„(u)d". 



TT 



vr 



n-l 



71-1 



r„(u)(i", 



n=l 



Theorem 15.11 and the definition of the norms yields 

an = 0{n-'), pn = 0(n-5), 7n = 0{n-^) and \\DE J^+i = 0((n - 1) Ivr"") . 
We apply to K29\l and obtain 

G = C-\Ei) + d^C-^{E2) + d^C-^Ei) + d^C-\E). (5.31) 
To the first two summands, theorem 14.41 applies and yields 

= 0((n-7)!7r-"), 



{C-\Ei)}„-a{n-l)\(^-J Tniu) 



n-l 



{C-\E2)}n-P{n-l)l 



vr 



Tn{u) 



\7r/ 



n-l 



0((n-5) Itt-""), 
0((n-3) Ivr-"), 



where 



^ OO ^ CXD ^ OO 

a = - "n, /? = - V 7 = - y] 7n 

vr vr vr ^-^ 



n=7 



n=5 



n=3 
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To the last part of (|5.29|) . we apply theorem 14.31 and obtain 



{C-\E)}n = 0((n-l)!7r-"log(n)). 



Using G = C (E), we have shown 



vr 



n-l 



G„-Q(n-1)! - Tn - /?(n-3)! - t„_2 - 7(n - 5) ! 



vr 



n— 3 



TT 



n— 5 



7"n-4 



o((n-7) Ivr-" log(n)). 



If we use (jS.Sp and the relation r„_|_i = y^DTn and we apply theorem 13.31 (5), we 
obtain our final result of this section 



i \ 71+1 
Fn-an\{—j t'^j^^ 



(7 \ n— i / 7 \ n— 3 

-) <_i-7(n-4)!(-) < 

o((n-5) Ivr"" log(n + 2)). 



and 



vr/ r2 



/3(n-l)!(^)""'^-7(n-3)!(- 
V vr / r2 V vr 

of(n-5) !^-"log(n + 2) 



T2 



6 Preparation of the functions to construct a quasi- 
solutions 

In the previous section , we have shown that equation (|2.4|) has a formal solu- 
tion and we found an asymptotic approximation of the coefficients of this formal 
solution. We will use this to construct the quasi-solution. To that purpose, we 
define the functions 

/ i\ n+l 

Hn{u):={n + l)\i^-j r„+2(u) (6.1) 



n odd 

We rewrite 



h{t,u):= Y,Hn{u)-—-. (6.2) 
n=7 ^ ' ' 



^(*'^) = (^) E (^ + l)(^ + 2)(^)"^-+3(n). (6.3) 

n=6 
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Using proposition (|3.ip -(4). we obtain 

Ht,u) = ^ 5: -(-) ^(^(tanh(O)), 

or equivalently 



n=6 
n even 



?i=6 ^ 



5(0 , 



where g{^) = — 2(tanh(0 — tanh(O^) and ^ = ^(u) = artanh{u). 
and hence that 



h{t, u) 



-1 



1 /it\n d 



vr^ z — ^ n ! V vr / d"? 

n=0 ^ 

n even 



2u 



where 

li{t,u) : 
Thus we can write 

h{t, u) 
This gives 

h{t,u) = (l-ti^) 



2 

3^6' 



-1 

2^L- 



it 



9{i + -) +9{i 



vr 



+ ^l{t,u). 



(6.4) 



(4n3 - 4u) sin d:)"^ - (Gu^ - 2u) sin (i)^ + 2u 



TT^f COS (-) + sin (-1 



(6.5) 



For fixed real u the function h{.,u) is analytic in |t| < p, where p = In the 
subsequent definition, we consider real values of n, < n < 1, here h{.,u) is also 
analytic with respect to t on the positive real axis. 
Now we define the function Ti.{d, u) by 



n{d, u) 



r+oo ^ 

:= / e~dh{t,u)dt, for (0 < n < 1). 
Jo 



(6.6) 



The function Ti.{d, .) is real analytic ; they can be continued analytically to the 
interval — 1 < w < 1 in the following way. Choose some positive number M and let 
Fi the path consisting of the segment from to Mi and of the ray 1 1— > t+Mi, t > 0. 
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Let T2 the symmetric path that could also be obtained using —M instead of M. 
Recalling (|6.4p . we can also define 



it 



r2 



vr 



it. 



e <ig[i + -)dt+ / e <ig{i--)dt 



IT 



where 



fii{d,u) :-- 



for — tanh(^M) < u < 1, where ^ = artanh(M), because the singularities of tanh 
are «(f + nir), n integer. As M is arbitrary, this defines the analytic continuation 
of n{d,.) for -1< n < 1. 

In the sequel we consider uq G] — 1, 0]. 

Lemma 6.1. // we consider the above function 7i{d,u) and the operators defined 
in (EH^- Then, foruo<u<l 

- 1. 

S2{n) = {i-u^)Vi{d,u), 

where the function n) is analytic, beginnings with d^. 

- 2. 

C2in) = -n{d,u) + {l-u^)V2id,u), 

where the function 'D2{d,u) is analytic, beginnings with d'^ . 

- 3 



1 



T2{u) 



n 



2ue 
T2{u) 



n{d,u) +V3{d,u) 



where the function T>3{d,u) is analytic beginnings with d^ 
Proof. (1)- If we replace u by r+ and T~ in (16. 5p we obtain 

n{d,T+) -n{d,T-) = J e—d(h{t,T+)-h{t,T')jdt, 
this implies 



n{d, r+) - n{d, t-) = -^{i+ - 1-) + a{d, t+) - cj{d, t- 



(6.7) 



where 



1- 
a{d, u) 



+00 



e-^(<7(e(T+) + ^)+5(C(T+)-5)dt, 



'e-2[g{aT-) + '^)+g{({T-)-'^))dt, 



+00 



e d ix{t,u)dt. 
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Using ^(T^) = ± d, where ^ and u are coupled by ^ = artanh(ti), we obtain 

r+°° t it /■+°° t it 
T+ = y e~dg(^^ + d+ -)dt + J e-dg[^ + d )dt. 

If we substitute t + nid in the first part, t — irid in the second part, we obtain 

. f+^+^^d , it. r+oc+nid 

J —nid ^ J -nid ^ 

We apply Cauchy's theorem 

T+ = -y^ e—d[g{i + -)+g{i--))dt + e—dg{^ + -)dt 

+ / ^"'9{C--)dt. 
Jo ^ 

Substitute t = —i s d in the second part, t = i sd in the third part, we obtain 
e--d(g{i+'^)+g{i-'^)yt + 2d sin{s)g{^ + ^)d s. 



This implies 

T+ = - r^e~Ug{^ + -)+g{^--))dt-Ad{l-n')V+{d,n), (6.8) 
Jo \ IT IT y 



where 



+ tanh(^))(l -tanh(^^2A 



V+(d,u)= sin(s)^ '-^^^ ' ^ ' U s. 

Jo {l + utanli{^)f 

We can also use the same method for I~ and obtain 

= -j^ e-d(g{i + -)+g{i--))dt-U{l-u')V-{d,u)A^-'^) 



where 



tanh(f ))(l -tanh(^)2) 



V'(d,u)= sin(s)^ -^^-^ ' ^ U s. 

Jo (i_„tanh(^))' 

Consequently 

n{d,T+)-n{d,T-) = ^^^^--^(v+{d,u)-V-{d,u)) +a{d,T+)-a{d,T-). 



Using (12221), we find 



S2in) = {l-u'^)Vi{d,u), (6.10) 
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where Vi{d,u) is analytic, beginnings with d^. 

(2) . Using the same method, we obtain 

n{d, r+) + Hid, T') = -2n{d, u) + 2(1 - u^)V2{d, u). 
With IKm we find 

C2{n) = -Hid, n) + (1 - u^)V2{d, u). 

T>2{d,u) is analytic, beginnings with d'^ . 

(3) . Using (12:22]) and (1), (2) of this lemma, we obtain 

M^'^) = '52(^)C2(W)+C2(^)52(H) 
\T2[U) J \T2\U)) \T2(u)) ^ ' 

2u£ 

= — --n{d,u) + V3{d,u) 

T2{U) 

where 

V3{d, u) := (1 + (1 + n^) smh{df) (1 - u'^)Vi{d, u) + 2ue{l - u^)V2{d, u). 

Proposition 6.2. We have 
1 . For uq <u <1 

OO -11 

~ + T„+2(u)d", asd\0, (6.11) 

n=7 ^ 
n odd 

2- \^{d.'^)\<Kd for UQ<u<l (d > 0). 
Proof. 



1. To prove (1) we use Watson's lemma and (16. 2p . if |i — /ctt^I < for uq < 
-u < 1 

2. - For Uq <u <\ 

■H{d,u)=—j e-*/'^h^-^\t,u)dt, 
where h^~^\t,u) satisfies 

|)V-) = |^ andf|)\(-)(0,n)=0, for = 0, 5. 



Because of 



h''-^\t,u) < K— fort>0. 
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we obtain 



This implies 



- for i < u < 1 



< Kd 



{d > 0). 



dH 
du 



(d,u] 



< K / e~ddt, 

Jo 

< Kd {d> 0). □ 



Let us consider a sequence Rn{u) of polynomials of degree at most n, such that 

\\Rn\\n = o(^{n-5) !7r-Mog(n)) 

Lemma 6.3. Let polynomial series difine R € Q with the above estimate be given 
and define 

r{t,u): = J2R^{u) . {teC,\t\<—,uo<u<l), 

n=7 ^ ' ' 



{t,u): = r(^,u) + {t-^)^{^,u), {t> ^,uo<u<l) 



TT 



2' dt"- 2 



n{d,u) : 



e d r{t, u)dt. 



Then 



1. r is continuously differentiable function on the set B of all {t,u) such that 

2 

u satisfies uq < u < 1 and t is a complex number and satisfies \t\ < ^ or 

2 2 

t > The restriction of r to uq < u < l,\t\ < ^ is twice continuously 

2 

differentiable. For fi/xed uq <u <1 the function r{t,u) is analytic in \t\ < 

2. TZ(d, u) is continuous, partially differentiable with respect to u, has continuous 
partial derivative and 



n{d, n) - ^ Rn{u)dP', asd\0. 

n=7 

3. \n{d,u)\ < Kd^ , \^{d,u)\ < Kd^ , for uq<u<\ (d > 0). 



(6.12) 
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The importance of our definition of TZ lies in a certain compatibility with 
insertion of the functions T^,T^ for u. First let 

oo oo 

n=7 n=7 

oo oo 

n=7 n=7 

We obtain a new sequences R^{u),R~{u) of polynomials of degree at most n. 
Theorem Ol and (Ol) imply 

\\R+{u)\l = o((n-5)!7r-"log(n)) 

\\RnHl = o((ri-5)!7r-"log(n)) 

Therefore we can use lemma [63] for R^{u), R~{u), and obtain functions Tl~^{d, n), 
7^-((i, n). 

Theorem 6.4. There is a positive constant K independent of d, u such that 

|7^+(d,n) -7^((^,^+)| < Kd^e'^d for (0 < d < do, no < u < 1) 
\n'{d,u) -n{d,T-)\ < Kd^e'^ for {0 < d < do,uo < u < 1). 



Proof The proof is exactly the one of [5]. 

Definition 6.5. Let T>[d, u) be a function defined for < d < do and uq < u < 1. 
We say that u) has property G if 

POO 

T>{d,u) = / e~^q{t,u)dt {0 < d < do, uq < u < 1) 
Jo 

is the Laplace transform of some function q{t, u) that has the following proprties : 

2 

L q{t,u) is defined if uq < u < 1 and either t is complex and \t\ < \- or t is 
real and t>Q, 

2 

2. q{t, u) is analytic in \t\ < ^ for uq < u < I, 

2 2 

3. q{t, u) restricted toQ<t<\ort>\-is continuous and the lim ^ q{t, u) 
exists for every uq < u < 1, 

4- there is a positive constant K such that 

\q{t,u)\ < Ke^^, for t > 0, uq < u < 1, {0 < d < do, uo < u < 1) 
The set of all function with the property G will be denote by Q 
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Lemma 6.6. For uq < u <1 

1. If7i{d,u) is the function defined in (6. 6]) then 

d^n{d, u) = {l- u^)n{d, u) + o(^{l- u^)e-^ 

where 7i{d, u) has property G 

2. Let k he a psitive integer. IfT>i,'D2 have property G and their first terms in 
the Taylor development at d = 0, begin with d'^ then 

Vi{d,u)V2{d,u) =d^V{d,u) + o[d^e~^ 

where T>{d, u) has property G 

3. If 7i{d, u) is the function defined in (6. 6|) then 



rL{d,uf = {I - u'^fd^£{d,u) + o{{l - u^fd^e'^^ 

where £{d,u) has property G 

4- Any function V{d,u) analytic in a neighborhood of d = has property G if 
V{0, u) = for all u, 

5. IfTZ(d,u) is defined by lemma \673\ then —^TZid,u) has property G 

6. //I?i,2?2 have property G then so do T>i + I?2> — ^2 and T>i ■ V2 

7. // ^{d, u) has property G then 

\V{d,u)\<Kd {0<d<^) (6.13) 
with some constant K > independent of u. 

Proof. 

1. For tt > 0, we have 

d^n = {l-u'^) e~'igi{t,u)dt (6.14) 
Jo 

where 



g-iit^u) = I I I I h{s,u)ds dr dv dO 



Jo Jo Jo 

g4{t,u) has a logarithmic singularity at tfc(s) = {2k + 1)— ± d^i for (A; > 

2 

0,s > 0). it is analytic in |t| < ^ and lim ^ gi{t,u) exists. 
If we put 



^{{dju) = / e~dg^{t,u)dt 
Jo 
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where 

f 2 

~ (4. \ J 54(^,^-'), if t < 



2 

then 7i(c?, u) has property G and 

d^n{d,u) = {l-u^)'H{d,u) + o({l-u'^)e-^^. 
(ii)- For — uo < u <1, where < < 1, we have 

'H{d,u) = / e~dh{t,u)dt + 2'!ri'^Res(e~dh{t,u),tk{s)] 

jfc>o 

TT vr 

where 2^ ^ ^- ^^'^ < ^ < '"o, this formula coincides with the formula 

/ e dh(t,u)dt 
Jo 

and extends it by real analytic continuation for —uq < u < 0. 
This implis 

d'^n{d,u) = (l-u^) J e-^gi{t,u)dt + o(^{l - u'^)d'^e-^^ 
we obtain 

d'^n{d,u) - {l-u^)'H{d,u) = {l-u^) j e-^g2{t,u)dT 
+ 0({l-u^)e-^), 

2 

where V is the path following the real line from infinity to (^,0), then along 

2 / 2 2 \ 

the vertical line from (^,0) to ( ^ tan(<y9) j and finally along the line 
y = tan(99).x from ^tan((/?)^ to infinity. 
Since g4{.,u) is bounded on F, then 

d'^H{d, u)-{l- u'^)n{d, u) < K{1 - u'^)e-^, 

where K is positive constant. Finally 

d'^n{d,u) = {l-u'^)H{d,u) + o(^{l-u^)e-^^ for (0 < < u < 1) 
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2. We assume that Vi{d,u),V2{d,u) have property G and their first terms in 
the Taylor development at d = begin with d*^. Then 

f°° t 
Vi= e~df{t,u)dt 
Jo 

POD 

V2= e~^git,u)dt 
Jo 

where f{t,u),g{t,u) are analytic in \t\ < ^ and f{t,u) = 0{t'^~^),g{t,u) = 

POD 

Vi{d,u)V2{d,u) = e-^f * g){t,u)dt (6.15) 
Jo 



Since 



h{t,u) = {f*g){t,u)=f f{t,u)g{t- s,u)ds 

Jo 



fit - s,u)g{s,u)ds 



t pt 
f{s,u)g{t - s,u)ds + / f{t - s,u)g{s,u)ds. 

t/2 Jt/2 



fit, u)giO, u) + /(O, u)git, u) - fii,u)gi^,u) 



2 

For t < the function hit, u) is k times differentiable with respect to t and 
h'it, u) = 

I fis,u)g'it- s,u)ds+ I fit- s,u)gis,u)ds 

Jt/2 Jt/2 

[ fis)g'it-s)ds+ [ f'it-s,u)gis,u)ds 

Jtl2 Jtl2 



+ 



1 1/2 

fiL'^)9iLu) 



h^''Kt,u) = f fis)g^^-'^^it-s)ds+ f f^^-^\t-s,u)gis,u)ds 

Jt/2 Jt/2 



n=0 



Observe that h^^\t,u) is continuous on [0, 7r^[, it is analytic for \t\ < \. If 
we put 



hit, u) 



hit,u), 

(t-^)\ ifi>4 
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then 

/ e~-dU^\t,u)dt, 
Jo 

has property G and 
iVi ■ V2)[d,u) = J e~dh{t,u)dt = J e~dh{t,u)dt + oi^d^e'^j 
= d!" j e—dU^\t,u)dt + 0[d^e-^ 
= d^V{d,u) + o(d''e-^ 



where 'D{d,u) has property G. 
3. We have 

n{d,uf = ^{d^n)-{d^n) 

Using (1) and (2) of this lemma, we obtain immediately the result. 

For < li < 1, the proof of (4), (5), (6) and (7) is exactly the one of [5]. This proof 
is valid for uq < u < 1. 

7 Approximate solution A{d, u) of the functional equa- 
tion ( 12:41 ) 



In this section we use the formal solution of (|2.4I) . The estimates of section 5 
for its coefficients and the preliminary results of section 6 allow to construct the 
function A{d,u) that almost satisfies (|2.4p . except an error which is exponentially 
small as d — > 0. 

In theorem 12.11 we found that ()2.4p has a uniquely determined formal power 
series solution 

A{d, u) = -B{d, u) = -iu+Y, ^n(n)<i" j , (7.1) 



where 



n=7 
nodd 



U{d,u) = eu + {u - u^)d^ + {^u^ - + 2u)d^ , 



(i \ n+l / i \ 

-) <+i(n)+/3(n-2)!(-j <_i(n) 

n-4 



+ 7(n-4)!(-) <_3(u)+i?„(n) 
= o((n-5)!7r-ilog(n)), (n > 7). 
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Consequently 



nodd nodd 

Now we set 

B(d, u) = (a + I3d^ + jd^) -^—-n(d, u) + TZid, u) + U(d, u), (d > 0,uo < u < 1), 

T2[U) 

(7.2) 

where TL^d, u) is defined in (|6.6|) and TZ[d, u) is the function corresponding to 
Rni n = 7, 9, .. By lemma [Ql and proposition 16.21 we have 

B{d,u) ~ B{d,u), as d\0 for every uq < u < 1. (7.3) 

Theorem 7.1. We have 

K TT^ 

A{d,T+) - A{d,T-) - f{e,A{d,u)) < — e'^, /or (0 < d < 4, lio < n < 1), 
where K is a constant independent of d, u and 

Aid,u) = ^B{d,u) (7.4) 



Proof. We set 

T{d,u) =S2{B)-e'^ + B{d,uf. (7.5) 

then 

T(d,u) = (a + pd^ + jd'^)S2(^rT^) + ^{d,u? 

\T2[U) J 

+ i7^(d,^+)-i7^(d,^-) + 52(^7) 

With (|7.2p and lemma [HTTl this implies 

T(d, u) = T{){d, u) ■ H{d, u) + Ti{d, u) + 7^((i, n), 



where 



ro((i,?x) = [a + (5d^ + -id'^) -2eu + 2U(d,i 
ri(d,K) = (a + (5d'^ + -fd^?^—:rn{d,uf + n{d,u) 

7-2('u) 

+ 2(a + /3d2 + 7d'^)^-7^((i,^i)7^((^,^i) 

T2{U) 

T2{d,u) = 2U{d,u)n{d,u) + ^n{d,T+) -^TZ{d,T~ 
%{d,u) = U{d,uf + S2iU)~e^ + ^V{d,u) 
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Observing that %){(!, u) is analytic function beginnings with . Using (1) of lemma 
?6l we obtain that dTQ{d,u)'H{d,u) has property G. 
Applying lemma ES, we obtain 



Ti{d,u) = d^V^{d,u) + o{d^e-^^ 



where V^{d^u) has property G. 

Using (5) of lemma 16.61 and theorem 16.41 we find 

T2{d,u) = d^Vi{d,u) + o[d^e-^ 

where Vi{d,u) has property G. 

On the other hand, —Ti{d,u) has property G. This implies 

T{d, u) = ^'^5{d, u) + K.{d, u), 
where V^[d,u) has property G and 

\K:{d,u)\<Kd^e-^ for (0 < d < do, no < n < 1). (7.6) 
Altogether we have proved that 

1 t 

T{d,u) = - e-dd{t,u)dt + lC{d,u), 
d Jo 

2 2 2 

where S{.,u) is analytic in |t| < it is continuous on [0, ^[ and ]^,oo[ and has 

2 

a Umit as t — > ^ for every uq < u < 1 and satisfaies 

\5{t,u)\<Ke''\ (7.7) 

with a constant K independent of u. Now if 6{t,u) = Yl'^=o ^n{t)t'^ is the power 
series of 6{t, u) near t = 0, Watson's lemma implies with (|7.6I1 that 

oo 

dT{d, u) ~ n ]5nd"'~^^, as d \ for every uq < u < 1. 

n=0 

On the other hand, because of its definition 

T{d,u) ~ ^B{d,T+)-^B{d,T-)-e'^ + B{d,uf 
~ + 0d + ---, 

because the formal series B satisfy (15.111 . But this means that all 6n must vanish 
and 5{t, u) = if 

|t| < ^, uo < u < 1. With (f7?7D . thus we obtain for uo < « < 1 

/"OO /"OO 2 

/ e-^\5{t,u)\dt < e^^Ke^^dt<Kde-^, {0 < d < do). 
Jo J^ 
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We obtain 

\T{d,u)\ < Ke~^ {0<d<do). 
Finally, we have shown that there is a constant K such that 

A{d, r+) - Aid, T-) - f{e, A{d, u)) 

for < d < do, < n < 1. 

8 Distance between points of manifolds 

Let y{t),t g] — l,oo[ an exact solution to the difference equation (12. ip . so that 
the stable manifold is parametrized by i i-^ [y{t),y{t + e)). We know that 
y~{t) = —y~{—t) is also an exact solution of (|2.ip with y~{t) — > 1 as i — > — oo 
. Morever the instable manifold W~ is parametrized by t i— > [y~{t),y~~{t + e)). 
We introduced the quasi-solution A{d,u) in previous section, then A~{d,u) = 
—A{d, —u) is a quasi-solution close to the exact solution y~{t). We denote by Wg, 
Wi the manifolds close to respectively W[~ parametrized by t i— > 
e)) respectively t i-^ + e)) , where = A[d,u{t)^ and = 

A-{d,u{t)). 

In this section we will show that the vertical distance between a point of the 
stable manifold and the manifold Ws is exponentially small as well as that between 
Wi and ; both are smaller, as we will show than the distance between Wg and 
Wi. This eventually proves that W^ and W~ do not coincide. 

In order to estimate the distance between some point {xQ^yo) on the stable 
manifold W^, we consider the sequence P„ = {xn-,yn), defined by Pn+i = (piPn), 
where 

Pn = {xn,yn) is a point of the stable manifold. We obtain 

\ynj \yn+lj \Xn + 2e{l-yi)J 

Since ^ = (pocp, where ^ is defined in (ll.3|l . and the stable manifold W^ is invariant 
by the application $, then the stable manifold W^ remains invariant by (p. Then 
theorem 17.11 implies 

^+{t + e)= ^+{t - e) + 2e(l - C+{tf) + \r,{t)e-^^ (8.3) 

where r]{t) is bounded. There is a sequence i„ g] — 1, +00 [ such that 

e+(tn -e) = xn= A(d, T- {d, uitn))) = ^(d, T~ {d, f/,(n))) (8.4) 



i 



T{d, u) I < — e 



2d 
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w 



(b , f(b )) 



Fig. 4 - The vertical distance between a point of the stable manifold and the approxi- 
mation defined using the quasi-solution A for the logistic equation . 



+ = >l(d,T+(d,C/,(n))) (8.5) 

where {/^(n) = u{tn). The relation (|8.4p implies that the point where 6„ = 

the vertical projection of the point on the curve t ^ 

e)). But, according to ([O) . (lO) . we have 

- e) = = Vn- (8.6) 

Now consider the error quantities 

de{n) = yn - A{d,Ue{n)) = yn - bn, (8.7) 

de(n) denotes the vertical distance between a point Pn and the manifold Wg. The 
manifold is parameterized by the curve 
{S,+(t),l+{t + e)). We obtain 

de{n + 1) = Vn+l - f{Xn+l) = Vn+l " fiVn) = Vn+l " f {de{n) + bn) , (8.8) 

where f{x) = ^+^^^^(x) + The point {yn,f{yn)) is the vertical projection of 

the point (a^ri+i, l/n+i) on the curve t i— > (^4-(t), ^+(t + e)) (See FIG. 4). 
If we use Taylor expainson, (|8.8p implies 

4(n + l)=2/„+i-/(6„)-/'(a„)4(n), < 0„ < 6„ + 4(n). (8.9) 

Then according to the definitions of / and 6„, we have f{bn) = C+itn + with 
(lO) and ([831), we obtain 

/(M = x„ + 2e(l - 6^) + ^r,{tn)e-^. (8.10) 

Together with jSJ]) and ([8?2]) this implies 
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ds{n + 1) = -f'{en)de{n) + 2e{bl - yl) - -??(t„)e-fe 



Using (18. 7p we obtain 



Now 



then 



d,{n + 1) = - [f'{en) + 2e{yn + 6n))4(n) - -r]{tn)e 



/(^0 = ii^f^, where = C+\9n) 



2e 



f\en) = 1 - 2u{Zn)d + [MZn? " + O(d^). 

On the other hand bn < On < Un imphes 

< 6„ < i{Zn) < Vn 

and thus 

f'ien) + 2e{yn + bn) > 1 + 2e{i+{zn) + bn) - 2u{zn)d 
+ (3u(z„)2-l)d2 + 0(d3). 

We know that 

i+{Zn) = u{Zn) + {u{Zn) - u{Znf)d^ + 0{d^), 

this with (I2.8|l imphes 

f'iOn) + 2e{yn + b„) > 1 + 2e6„ + 0{d^). 
As 6„ > 0, then for sufficiently small d, 

fiOn) + 2e{yn + bn) > 1. 

This implies 

f'{dn) + '^eiVn + bn) = 1 + 96 {u{zn)) , with some ge {u{zn)) > £C, 
ge{u{zn)) = 0{e) and c is a positive constant. With ()8.12|) . this implies 

de{n+l) = -[l + ge{u{zn))jde{n) - -r?(t„)e~2F. 

As de{n) ^ as n ^ +oo, this implies that 

l*(«)l < \e-i ± 

k=n 



ur'(^+9siu{z,) 
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provided this series converges. This follows from 



E 

k=n 



Visk) 



+00 



k=0 



where we used that rj is bounded. This imply that there is a constant K > such 
that 



k=n 



ViSk) 



< 



K 



With (183|) imphes 



and in particular that the vertical distance (ie(0) of (xo,yo) to Ws is 

4(0) = o(lexp(-g). 



(8.16) 



Now we will calculate the vertical distance between the twomanifolds Wg and 
Wi. To do this, we will need to extend both quasi-solutions in an adequate way. 
We define 



i+{t) = A{d,u{t)), for -^<n<l, 
l_{t) = A'{d,u{t)), for -l<u<\^ 



(8.17) 



First, we evaluate the quantities 

de{i+{t),L{t) 
using JTJ) and (TTiD 



1 1 

A{d,u) - A~{d,u), for - - < n < - (8.18) 



|+(t)-|_(t), (8.19) 
- (a + (5d^ + 7d^) {n{d, u) - H' {d, u)) , 



where 



Hid, u) 

n-{d, u) 



_s 11 

e <ih{s,u)ds, for — - < u < -, 



e d,h(s,—u)ds, for 



1 1 

- <u< -. 

2 - - 2 



(8.20) 
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Modifying the path of integral, we obtain 

i-L{d,u) = j e~dh{s,u)ds + 2TTiRes(^e~dh{s,u),Sk{t)^ 

^ lm{sk{u))<0 

— 2TTiRes(^e~dh{s,u),si.{t)^, 

lm{sif:{u))>0 

where Sk{t) = ^ ± ^ i + kir^ for A; > 0, F is the sum of two paths Fi and F2. 

Here Fi consist of two segments of the higher half-plane, one of these two segments 

2 

is parallel to the axis y = and begins at the point 1), the other ends in the 
point (0,0), and the path F2 is the symmetry of Fi relatively to the axis y = 0. 
This implies 



~ , , , /■_£,, s , TTc"^ ' ^ ( (A;-|-l)7r 

H{d,u) = J^e dh{s,u)ds + 2^2(1 _^2) Z^^^P ~ 



2d 



k=0 

ire-^' ^ f {k + l)7r^ 
+ 2^2(1 -u2) 2^ ^^Pl^ 2d 

Therefore 

n{d,u) = e-ihis, u) ds + g exp ( - ^^^) ■ (8.21) 

Similary for 7i.~{d,u) 



f s 7rcos(— ) 

'H~{d,u)= / e~d/i(s,ii)(is - — ^y^exp 



2d 



fe=0 

Therefore 

,/;r/x;r/x\ 27r(a + /3d2+7d4)cos(^) / vr^X ,^/cos(^) a.^X 

for — arctanh(^) <t< arctanh(^). Consequently 

/- , - , A 27racos(— ) ^2 /I ,7r^ \ 

Now we want to estimate the vertical distance between a point on the quasi- 
manifold Ws and its vertical projection on the quasimanifold Wi. For this purpose 
let t such that — ^ < u{t) < ^. Then the point ^^+(t), |^_|_(t + e)^ is on the qua- 
simanifold Ws- We suppose that (^^+{t),^~ {ti + s)j , where — | < tt(ii) < |, is 
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its vertical projection on the quasimanifold Wi, then the vertical distance between 
these two points is 

dist^it) = i+{t + e)- i-ih + e)= i+{t + e) - (8.23) 



where f{x) = i-[iz^{x) + e j , the inverse function of ^_ . Using (|8.22p and 
we obtain 

dist^it) = i+{t + e) - f{i^{t) + ee(t)), 
where ee(t) = ( ) . We use Taylor expansion and obtain 



dist,{t) = + e) - f{L{t)) - ee{t) ■ f'{9), 

= i^{t + e)-t{t + e)-es{t)-f'{d), (8.24) 

where |_(t) < 6* < |_(t) + Ceit). With IKOB this implies 

dist^it) = eS + e)-ee{t)- ~f'{e). (8.25) 



, X . t(t-\0)+e] 
Since Ceit + e) = -eJt) + Oie ■ eJt)] and /'(0) = — ^ ^ = l + 0(e), 



we obtain 

disUit) = {-2 + 0{e))ee{t). (8.26) 

Consequently 



disUt)r^-— —-^^6-27. (8.27) 



47racos(^) 
e3(l-tanh(fi)^) 
Combining (I8.16P and (|8.27p . we conclude 



,2 



9 An estimate for the main asymptotic coefficients 

In this section we show that the number a in the main theorem of the introduc- 
tion is not zero and we want to find bounds for an and will then obtain estimats 
for a using the first a„ which can be computed explicitly. We recall that 



00 

i 

a 



^ 00 

(9.1) 



71" 

n=7 
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where an{n — ^n" is the leading term of the coefficient {C{G)}n{u) of 

in C{G){d, u), where G = Q{d, u).F{d, u), F{d, u) = B{d, u) - U{d, u) (cf.K^ 
and KB)- 



For a polynomial series X{d,u) = Yl'^=o ■^niu)d'^ where the degrees of Xn(u) 
do not exceed n, we can write Xn{u) = ^^=QXnkU^- Then we denote by X{z) = 

The mapping X{d,u) i— > X{z) extracts the leading terms of the series X{d,u). 
It is compatible with addition and multiplication. If we define the following ope- 
rator D = —z'^-^, we also find 



DX{z) = DX{z) = -z^X'{z) 
C{X){z) = cosl.(^\x{z):=\(x(^\+x{ 



2 1 ' ' 2\ \2 + zJ \2- z 



C^){z) = cosHD)Xiz):=Ux(-^)+x' ' 



2\ \l + zJ \l-z 
1/;^/ z \ '-z z 



S^){z) = sinh(5)X(z):=-(l(^)-X^^_^ 



where the operators 5,C,52,C2 are defined in (I2.18p . If we apply the mapping 
X{d,u) I— > X{z) to the equation (|5.13l) , we obtain 



c(^I-yz)=n{E){z), (9.2) 
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where 



C{x) 

n{x) 



E 

y 

Q 
Qi 

Vi 



W2 



Wi 



'^\y + ^c(^]y^ 



QiJ 2z^Q 222Q2 



W2^fy\ 



= C(G)= £ a„(n-l)!(- 



n=Oodd 



n-1 



-z' + z'- 
-z" - lz\ 



-z^ + — 
6^ +18^ ' 



- 192 - 240^2 + 508/ + 1995^^ + 3418^^ 



+ 4389z^° + 217z 



a2 



3(^2-4)4 



.11 



591 - 15457^2 + 21610z^ + 59194^^ + 36583z^ 



+ 2759z^° 



27(^2 - 1) 



^11 



2385 - 32832^2 - 11790^^ + 231917/ - 116267/ + 720682^^° 



+ 1018936^^^ + 290089^^^ + 2480z^^ 
Wi = VC2{Q) + WMQ) 



81(/ - 1)^ 



= —z 



54 - 432z^ + 1449z^ - 2724z^ + 4449z^ - 14719^^° - 12857^^^ 



+ 951462^'' + 110905z^^ + 29041z^^ + 248z 



W3 + PV3, 
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54(/ - 1)^ 



= -8/ 



- 1152 + 1168/ + 9744/ + 18901/ + 14664/ 



+ 1240z^° 



3(/ -4)^(4 + /) 
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5270002^° - 1988300z^^ + 14763570z^^ - 254881452^'' 



+ 11808105z^2 + 113810272^° - 17065562z^^ - 15479497z^^ - 15777878^^'' 
- 21499321z^2 + 14844222z^° + 4847178z^ - 1419003^^ - 5410262"^ 

972(^2 - 1)^2 , (9.3) 



+ 320490z^ - 50220 



with X = x{z) is a power series. 

Now for integers m > p > —1, m > 0, we introduce the set 

J^{m,p) = < ^^*x„z"'\ x„/ [(n — p) Ivr""] , n = m, m + 2... is bounded > 

^ n=m ' 

M.{jn^p) is a Banach space with the norm 

|a;|m,p = sup { [(n — p) Ivr"*^] , n ^ m, n = m mod 2}. (9.4) 

oo 

Here the symbol for the summation means that only summands with 

n=ra 

n = m mod 2 are considered. 

Proposition 9.1. Let m,p,q be integers and m > p > —l,m > q > — l,m > 0. 
The norms ^9.4\ ) have the following property : 

1. If q < p < m then M{m,p) C M.{m, q) and 

\x\m,q < 7 ttI^U-p /^^ ^ ^ M{m,p) (9.5) 

[m — q)\ 

2. If X e M{m,p) then zx £ M{m + l,p + 1) and 

\zx\m+l,p+l = TT\x\m,p (9.6) 

3. D : M{m,p) i-^ M{m + l,p) and 

'n-\x\m,p < \D{x)\m+i,p < vr —\x\m,p (9.7) 

m — p + 1 

Now we consider an even power series f{z) = X^^o*/?'-^"' having a radius of 
convergence greater than vr. Then f{D) maps ^A{m,p) into itself and 

\fiD)x\m,p < \\f\\m,p\x\m,p, foix £M{m,p) (9.8) 



where 



sups ^__^* l/il'^"' ]__[ ■ — -\n>m,n = m mod2>ifp>l, 

j=o k=l ^ J 

sup < l/jlvr-'^n > m, n = m mod2>ifj»<0. (9.9) 

I 7=0 J 



m,p 



m,p 
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We have also 



\f-x\m,p < Mm,p{f)\x\m,p, for X G M{m,p), 



(9.10) 



where 

MmAf) = sup I \^ >m,n = m mod 2| . (9.11) 

We conclude that the inverse of C maps Ai{m,p) with p > 3 into ^A{m, 1) and 

(9.12) 



|C (3^)|m,l ^ Cinp\x\mpj With Cmp — / 

vr — ' 

n=m 



{n-l)V 



Let X € ^A{m,p) and y € Al(m, q) and assume that m + g < n + p. Then we 
obtain x.y G ^A{n + m, q + m) and 



(9.13) 



with 



7V- 



ar,s = sup 



>r + s,N — r — s even 



Finally we need estimates for the inverse of the linear operator C in (|9.2p . £ maps 
A^(m — l,q) onto M.{m,q) and is one to one. Now let x G 7W(m — l,p) and 
w = C{x). If we proceed similar to the proof of theorem 14.61 we find that 



Dy = f{D) 



(4+_£^ 



-w 



(9.14) 



where f{z) = zj sinh(z/2). 

We can apply (|9.6p .. (19.131) and obtain for m > max(p + 1,4) 

iw)\m-l,p < Lm,p\w\m,p foTW G M{m,p), 

where 



Lm,,p 



1 



vr 



sinh(z/2) 



m— 2,p— 2 



(4 + z^) z^ 
4^1 



(9.15) 



(9.16) 



Now we are in a position to use (19.21) to estimate E. We use that E can be 
computed recursively from (|9.2I) . If a; is an odd power series and x = E mod 
with n > 9 then via C~^(x) = G mod z"- we first obtain TZ{x) = 1Z{E) mod 2;"+^ 
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and then L ^ [TZ^xyj = -M- mod z^+^.Thus if xq = mod is arbitrary and we 
let 



Xn+i ■■= QiC ^(lZ{xnfj forn = 0... 



.17) 



then Xn = E mod z^^^" for all n. We begin with the polynomial xq of degree 
< 79 with Xq = E mod z^^ . We computed it explicitly and found that 



\x0\9 7<Ko:= 308027.359777894414. 



(9.18) 



Then we consider the sequence x„ of polynomial series defined in (|9.17|) . which 
now satisfies Xn ^ E mod 2;8i+2n-_ -yVe will show 



a^nlgj < for alln > 



(9.19) 



and this implies \E\<^j < Kq because each coefficient of E is also that of some Xn 
if n is sufficiently large. 

We prove (I9.19P by induction. For n = it is true. Assume that [xnjgj < Kq, 
for this Let us take again the equation (I9.18P with (|9.2p . We obtain 



Xn+i + Q\C: 



1 ^ 



-1 i W^n( ^ 

Qi 



+ 



\2z^Q 



\Q 



(9.20) 



where = C'^ixn). Then using (lOl) . . (l9T7ll . we obtain 



\2z^Q , 



2z^Q^ \Q 



11,10 



17,10 



17,10 



15,10 



< ^lko|9,7, 

< ^2ko|9,7, 

< K^\xq\<)j, 

< K^lxollj 
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where 
Ki : 

K2 : 



cosh(|) 



7,5 



^')||cosh(z)||^_^M9,i^4- ) Cgj, 



2z 



4 • 



Ml3« 1 + -z' Li4,8Mi4,8l TTT 1^18,12 



W4 

2^ 



Q 



2 ,"8,8 C9_7. 



This imply 



Xn+l + QlC ^ ( ^ 



11,10 



< (E:i+i^2 + ^3+i^4i^0)i^0 



Evaluation of these quantities yields approximately 





K2 






375142.8501573 


633.1622891589 


10214.44411459 


0.000011150281 



and altogether 



Xn+l +QlC ^{ ^ 



< 118896679183. 



(9.21) 



11,10 



Now we consider g{z) := z C y^j , letting z y = z g{z), with (|9.15|) we 
obtain 

-/}(z)(4 + z2) 



where 



Dy = D{z-^g)= f{D) 



Az^Vi 



(9.22) 



sinh(z/2) 



The coefficient of z^ in f{z) is smaller than 4(27r)^ we obtain by comparison 
of the coefficients of z^ in (|9.22p 



0--l)|g.+3l<4i:*(2vr)'=-^/i.|-^ 
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and hence 

ipl — ^ < 4 V -p — '— < 3668333. 
0-2)! - ^ {k-l)\ - 

We obtain for add j > 9 



(j - 10) 

and 



. < vr^ . 3668333 • (j - 9)(j - 8)(j - 7)(j - 6)(j - 5)2 



-i+3 



ij-2\ 



(i - 10) iTT-J - ■ ^^^^^^^ ■ " ^^^^ ~ ^^^^ ~ 

but, |^g(-2^)| . = gj + igi-2i with (|9.20|) this yields eventually 

\{^n+l}j\ 



(j - 10) In-i 

This implies 



< 118896679184 for j > 81. 



(j - 7) Ivr-J - 72 • 73 • 74 

and with Xn+i = Xn mod we finally proved that |x„+i|9j < Kq = 308027.35. 
Thus the proof of (|9.18l) is complete. 

Theorem 9.2. 

98048.15 

- {n- l)(n - 2)(n - 3)(n - 4)(n - 5)(n - 6) J^^"^-^' 



and even better 



(jl — IQ) I 

\a„\ < 37845988419-^- ^ for odd n > 81. 

(n — 1) ! 



As a consequence we obtain estimates for a 

N 

(n - 1) ! 



1 ^ 



a — 

n=l 



< 12046752267 V ^ — for n > 9, 



n=N 



< 12046752267 • ( 1 + 



< 



16(iV-2)! V N-1 
762333542 



A^(A^ - 1){N - 2){N - 3){N - 4){N - 5){N - 6){N - 7) 
We evaluated for n = 9, 11, 81 and obtain 

\a - 1.264150331[ < 6 • 10"^. (9.23) 
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